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Abstract 

A first-order perturbation approach to fe = Friedmann cosmologies 
filled with dust and radiation is developed. Adopting the coordinate gauge 
comoving with the perturbed matter, and neglecting the vorticity of the 
radiation, a pair of coupled equations is obtained for the trace h of the 
metric perturbations and for the velocity potential v. A power series 
solution with upwards cutoff exists such that the leading terms for large 
values of the dimensionless time ^ agree with the relatively growing terms 
of the dust solution of Sachs and Wolfe. 



1 Introduction 



The classic prediction for the temperature fluctuations of the cosmic background 
radiation by Sachs and Wolfe |^ overestimates the experimental value Q by at 
least two orders of magnitude. Worse than that, this prediction has been ob- 
tained by neglecting the temperature fluctuations on the surface of the last 
scattering, and including only the gravitational perturbations along the subse- 
quent path of the photon. Obviously, if the initial fluctuations are random, their 
inclusion can only increase the effect. While these computations can be criticised 
on the basis that the results are not invariant with respect to the choice of the 
initial hypersurface where the photons originate this is not likely to be the 
way of improving the results. Meanwhile, the discovery of large-scale structures 
(such as voids or walls) certainly did not bring us closer to the solution of this 
cosmic puzzle. It is hard to take seriously the suggestion[P that the observed 
fluctuations are primordial (as opposed to propagation effects), unless we are 
able to reduce the magnitude of the Sachs- Wolfe estimate accordingly. The way 
the surface of photon emission is deflned is just one detail of the gauge choice. 
As we seek diminuation of the effect, we may as well stick to the comoving 
gauge. For gauge invariant treatments, cf. Refs. 

Sachs and Wolfe have obtained, in a closed form, the first-order perturba- 
tions of a Friedmann universe with a flat 3-space, and filled either with dust or 
radiation. They assumed that the domain of the Universe in which the pho- 
ton travels is matter-dominated, and computed the temperature fiuctuations in 
the dust-filled universe. Their assumption is justifiable because the decoupling 
occurs near the time of equal matter-and-radiation density. (Some authors es- 
timate that the two phenomena occur simultaneously!^ while others take 
that the equal-density epoch precedes the surface of last scattering.) According 
to the unperturbed models, the radiation density pr dies out faster than the 
matter density, pm ■ However, due to the instabilities, the dust dominance may 
not hold everywhere in the perturbed models. 

Here we consider the refinement of the Sachs- Wolfe scheme by computing 
the perturbations of the fc = Friedmann universe in the presence of both 
dust and radiation. Instabilities are known to exist in two-fiuid cosmologies Q. 
Assuming, for instance, that the observed large structures are lately emerging 
manifestations of the instabilities, one might be able to reduce the magnitude 
of the metric fluctuations affecting the photon orbits. The energy-momentum 
tensor is a sum of those of the two media, 

^ b ^ ^rnb + ^rb K'-j 

The contribution of the dust has the form 

r^, = -p„,w'^wfa . (2) 
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For the radiation, 

T;'b=-lprU''Ub + ^PrS1 (3) 

The four- velocities are normahzed u^Ua = 1- We adopt the conformal form of 
the metric 

gab ^ {V) iVab + hab) , (4) 

where the unperturbed {hab ~ 0) metric satisfies |pl| 

3l('^V-^^-^=0. (5) 
\dri J a 

such that the density pi where i stands either for m (matter) or r (radiation), 
equals pmo Pmr ^-t somc prescribed conformal time t] = i^q. The solution has 
the form 

a = ^Xt]'^ - p (6) 

where the constants are defined 

\ ^ 3 Pro /-V 

2 The perturbed model 

In the perturbed universe, hab 7^ 0, the densities of the components can be 
written to first order 

= + Sp, . (8) 
Here pr and pm are the unperturbed densities 



Pniiv) = PrnO^ Pr{v) = PrO^ (9) 

and Spi are the first-order density perturbations. The indices of the perturbed 
quantities are lowered and raised by the Minkowski metric r/ab = ?7"^- 

It can be proven]!^ that a homogeneous universe cannot develop perturba- 
tions with comoving dust and radiation. In fact, the dipole effect of the cosmic 
background radiation provides an experimental value for the local relative ve- 
locity of the order of 100 km/sec. Thus we have good reason to assume that 
6ul^ ^ Sul- We then choose coordinates comoving with the matter: 

< = K (10) 

< = u^ + K (11) 
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where the coincident unperturbed velocities are 



< = ^'^0 ■ (12) 

The normaUzation conditions imply that /iqo = and Su'^ — SurQ — 0. 

After the decoupling, we may assume that the conservation laws apply sep- 
arately both to the matter and the radiation components: T^^.^ = and 
-^rb a ~ Thus Coupling occurs only via the universal gravitational interac- 
tion. We get two energy conservation equations from the timelike components 
and two momentum conservations laws from the spacelike components. 

2.1 Energy conservation 

The first-order perturbations of the dust satisfy the equation 

Pm ^ J 

where /i = ft,^ is the trace of the spacelike perturbation of the metric and a prime 
denotes partial derivative with respect to the conformal time 77. The solution is 

[E{x^)-h^ (14) 

where the integration function E{x^) depends only on the space coordinates 
a;"(a =1,2 or 3). 

For the radiation, the energy conservation law has the form 

(plV5<)^ = (15) 

or 

'J^ + ^/^)' + ««),.=0 (16) 

2.2 Momentum conservation 

The conservation law for the dust has the form 

{ahoj = (17) 

with the solution 

ah^o ^ F^ix") . (18) 
The remaining coordinate freedom makes it possible to arrange [|l3[ 

= . (19) 
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For radiation, the first-order momentum conservation law reads 

[pr'^'SUr^y ^\ap-'^Hpr^^. (20) 

The perturbed Einstein tensor wiU be written 

G% - oG% + 6G% . (21) 

Here oG% is the unperturbed and 6G% the first-order part. The Einstein equa- 
tions for the first-order quantities are 

6G\ = -{5pr + 6pm) (22) 
(5G"o = -^aprK (23) 



= liS^Spr (24) 



Substitution of the Sachs- Wolfe |^ expressions for 5G\ and separating the trace- 
free part of the metric perturbation 

Sap = - ^?7q/3^ (25) 

yields 

S^'\^, + ^A/i - 2^/i' = - 2a2 (4,, + 4,^) (26) 

S''\'- Ih'"' + A/jO" _4|2^-— -^ap.< , (27) 

6 \ a j A 

2h" + 4^h' - S^"^^, - ^A/i ^ ~2a'dpr . (28) 



.a' 
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-^^"^-^'^"^A/i (29) 
Taking the sum of (|2^) and ( p8|) we get the simple relation 

/i" + ^/i' = -a2(2(5p^-H(5p™) . (30) 
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3 The velocity potential 

For economy of writing, we introduce the scaled radiation velocity perturbation 

Va{Xa,ri) = pl^'^SUra ■ (31) 

(Note that the component vq = 0.) From the 77 derivative of ( [l^ ) and eliminating 
the term {^r) "-"^ ^he radiation momentum conservation ( ^o|) we get 

the relation 



/I ' 



(32) 



Taking the time derivative and the gradient, and using again ( pOD to get rid of 
the divergence term, the left-hand side is still a total time derivative. Integration 
and comparison with (|3^) leads to the uncoupled equation 



(33) 



where the integration function Qf}{x) depends only on space coordinates. We 
decompose the three- vector Va 



where the vorticity cua is divergenceless 



(34) 



(35) 



and V is the velocity potential. The terms containing the velocity potential v 
cancel in Eq. (|3^). Thus we get the simple relation for the vorticity 

Aup = -Q^ix) . (36) 

Following a suggestion of , we shall henceforth take the trivial solution uja = 
and Qp = such that we have 

Using the velocity potential in Eq. (|2^), we get 



(37) 



(38) 



This can be integrated. Since only the gradient of the velocity potential has 
physical meaning, the integration function V{r\) may be chosen to vanish. Thus 



1 



rap, 



(39) 
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Hence, Eq. (Hq) can be rewritten 



3w" -Av + ^pli'^aoh' = . (40) 



Introducing the velocity potential ( p9[ ) in (|30|), we get 

h" + -h' + p„,o- ( E{x) - h] + 8apl/\'' = (41) 
a a \ 2 J 

We may get rid of the inhoniogeneous term E{x) by introducing the function 

f^h- 2E{x) . (42) 

The explicit form of the scale factor a can be made simpler by use of the di- 
mensionless time variable 

f = i/^,. (43) 

Then we have 

a = n(e-l) f^^^ao) (44) 

such that the Big Bang occurs at ^ = 1. Eqs. ( ^ ) and ( pl[ ) take the form 

v-LAv + ^f ^ (45) 
K 

{e-l)f + 2U-&f + K-^v = 0, (46) 

where an overdot denotes derivation with respect to the dimensionless time ^ 
and the constants K and L are defined by 

K = 16V3ft„oA:„'''' (47) 

These coupled linear equations for / and v have an elaborate structure, even 
in the asymptotic regime. In the special case of the velocity potential being 



stationary, u = 0, Eq. (46) becomes the n — 2 Legendre equation for /. In the 



generic case, however, the expansion in Legendre series is cumbersome. 

4 Recursion relations 

The power series expansion in the time variable ^, 

oo 
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v{x,0 = E ^"(^)^" (50) 

n— — oo 

with Eqs. (45) and (EQ), yields the recursion relations for the coefficients: 



f6 

{n + l){n + 2)bn+2-LAbn + —in + l)a„+i = (51) 

K 

{n + l){n ~ 4)a„_2 - (2n^ - 6)a„ 
+ (n+ l)(n + 2)a„+2 + /s:(n+ l)6„+i = (52) 



Eliminating a„ we get 



{n + l)(n - l)(n - 4)6„_i - (n + l)(2n2 + 10)6„+i 
+ (?! + l)(n + 2){n + 3)6„+3 ^(n + l)(n - 4)Afe„_3 

71 — 2 

^A6„_i -L(n + l)A6„+i = 0. (53) 

n 

The convergence of the series expansion deserves further investigation. The 
terms with negative exponents will decay with time, and it is expected that 
their contribution to the temperature fluctuations becomes negligible. Hence, 
to evade fast growing terms, it appears worth while to seek solutions whith a 
maximal value of the exponent n. 

It is clear from the structure of Eqs. (^ij) and (^) that a cutoff of the power 
series can occur only simultaneously in a„ and 6„. Further necessary conditions 
for cutoff may be obtained from Eq. (|5^). The highest possible nonvanishing 
coefficient is bi. From ( [5l| ) it then follows that 02 is the highest coefficient of /. 
It is impossible to have a cutoff simultaneously both upwards and downwards. 
Our choice of the velocity potential implies that any coefficient 6„ satisfying 
A6„ — must vanish. Thus the leading terms of the solution of Eqs. ( |5l| ) and 
TTl) are 



+a_3(x)^ + ia_2(a;)^ + ... (54) 
v{x,^) = bi[x)i + b-:i{x)^+b-i{x)^ + ... (55) 
where the coefficients a„ and 6„ satisfy the relations 

ao(x) = -^a2 Abi{x) ^ ^a2{x) 
a_2(x) = fb,{x) Ab^six) = f 61 
a-4(x) - ^ Ab^^ix) = ^a.^ix) . 
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The leading terms in the density contrasts are 

^ = Ap-J'^a^^ {b,{x) 3&_3(x)^ - 46_4(:r)l + (56) 

-i ^02(2;)^^ +ao(a;) +a_2(a;)^ +a_3(a;)^ + ...^ ■ (57) 

Thus we establish the comforting result that for large values of ^, the leading 
power in time of the metric function h coincides with the power of the relatively 
growing mode of the pure dust solution [Q. For a more detailed account of our 
approach, c/. ]T5[ |. 
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